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TECHNICAL MEMORANDUM 


AN IN-DEPTH PROBABILISTIC STUDY OF EXTERNAL TANK ATTACH RING 

I. INTRODUCTION 


This work is part of an overall study in the area of evaluating the performance of the shuttle. 
Besides the orbiter and space shuttle main engines (SSME’s), the other main items of a shuttle are the 
solid rocket boosters (SRB’s), external tank (ET), and the ET attachment (ETA) ring. The present study 
is restricted to the ETA, but the concepts discussed are general in nature and can be applied to any 
structural component. The strengths of steel, the cross-sectional dimensions of the steel section, and thus 
the strength of the model itself and the loads are all variable parameters. In other words, taking the entire 
structural element under consideration, the resistance and the applied loads are random variables. 
Consequently, the traditionally used deterministic analysis approach (the so-called factor of safety) does 
not reveal the actual safety reserve in the structure, because it does not include the inherent variability in 
the material, geometric, and load parameters. Such a variability can easily be included in the probabilis- 
tic analyses using probabilistic methods which have been in use and continually improved over the last 
two decades. The objective of this research work is to use some of the existing probabilistic methods to 
calculate the reliability of the ETA ring at various critical sections for the limit state of stress. This is 
done both in terms of the traditional probability of failure (/?/) and reliability levels as well as the well- 
known safety indices (/J) which have become a commonly accepted measure of safety. 

For any probabilistic analysis, a certain deterministic model is needed for which the limit state 
function can be formulated in terms of the general basic variables of resistance and load. The variability 
of the materials, cross-sectional properties, and load can then be incorporated into the deterministic 
model, and the p value can be calculated which forms the basis of the probabilistic model. The detail 
work of the deterministic model and the corresponding deterministic analysis were done by United 
Space Boosters, Inc., (USBI) 1 under a general contract to Marshall Space Right Center (MSFC) and 
form a basis for the present probabilistic study. 


II. DETERMINISTIC MODEL 


The ETA ring used in the present analyses is part of the solid rocket aft booster assembly and is 
located at SRB station 15 1 1.0. It is the interface between the ET attach struts and the SRM. The ETA 
ring is shown in figure 1. It extends a full 360° around the SRM case, and its cross section is shown in 
figure 2. The details of the ETA ring are given elsewhere. 2 The pertinent details needed for the prob- 
abilistic study (discussed later in this report) are briefly dealt with in this section. The primary purpose 
of the ETA ring is the distribution of concentrated attach strut loads circumferentially into the SRM 
case. The strut attachment is a pinned connection designed to react loads in the plane of the ring. 
Secondary out-of-plane loading exists if the struts are misaligned. Four splices are required for the ETA 
ring configuration in figure 1. A typical spliced section is shown in figure 3. The outer systems tunnel 
splice plate is 4340 steel, heat treated to 180 lb/in 2 minimum tensile strength. The inner splice plate is 
4130 steel, heat treated to 180 lb/in 2 minimum strength. 

The critical sections of the deterministic model of the ETA ring are found to be the tunnel splice 
and H-fitting. 1 These are shown in figures 4 and 5. Figure 4 shows the strut fitting, and figure 5 shows 
the tunnel splice. 




Figure 1. SRB/ET attach ring (insulation omitted for clarity) (from ref. 1). 







vD 

Z 

P 


MOTORXASt 


Figure 3. ETA ring splice section at 154°, 289°, and 342° (from ref. 1). 
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Figure 5. Systems tunnel splice (from ref. 1). 


in. DETERMINISTIC FINITE ELEMENT ANALYSIS 


Finite element analysis has been used 1 to analyze the ETA ring using the powerful finite element 
code NASTRAN. The finite element model is shown in figure 6 which actually shows the underformed 
plot of the ETA ring. It was developed using MSFC/NASTRAN, version 65A on the TBE VAX 1 1/780 
computer system and executed on the Cray system at MSFC. 


Z 



Figure 6. ETA ring NASTRAN model boundary conditions (from ref. 1). 

The ETA ring finite element model is described in detail in references 1 and 2. The model con- 
sisted of approximately 10,000 nodes and 20,000 elements including some of the elements which are 
quadrilateral and bar elements. The ETA ring modeling details are shown in figure 7. The boundary 
conditions imposed on the NASTRAN model are shown in figure 6. NASTRAN output was obtained in 
the form of a punch file for the 18 “unit” load cases described below. 

1. SRM radial internal pressure (100 lb/in 2 ) 

2. Upper strut aligned (1 kip tension) 

3. Upper strut aligned (1 kip compression) 

4. Upper strut misaligned 90° forward (1 kip) 

5. Upper strut misaligned 90° aft ( 1 kip) 
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Figure 7. ETA ring modeling details (from ref. 1). 
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6. Lower strut aligned (1 kip tension) 

7. Lower strut aligned ( 1 kip compression) 

8. Lower strut misaligned 90° forward (1 kip) 

9. Lower strut misaligned 90° aft (1 kip) 

10. Diagonal strut aligned (1 kip tension) 

1 1. Diagonal strut aligned ( 1 kip compression) 

12. Diagonal strut misaligned 90° forward (1 kip) 

13. Diagonal strut misaligned 90° aft (1 kip) 

14. SRM axial load (10 6 lb) 

15. SRM bending moment M y ( 10 6 in-lb as sinusoidal axial load) 

16. SRM bending moment M z ( 10 6 in-lb as sinusoidal axial load) 

17. Web/IEA cover axial pressure (1 lb/in 2 ) 

18. IE A cover radial pressure (1 lb/in 2 ). 

FORTRAN postprocessing codes were written to multiply the “unit” results by the actual loads 
and linearly combine the results to create stresses and forces for each element and each load case. As 
pointed out earlier, the critical sections to be studied in the probabilistic analyses (based on the determin- 
istic finite element analyses) are tunnel splices and H-fitting. The load conditions studied are: 

USLI108A, USLI045A, USHI005D, and USHI006D. For example, for USLI108A, the corresponding 
equation for the total stress or force R T (for 0° alignment) is given by: 

Rj= 8.55/?i+9.15/?4+0.0/?i7+0.0/? 1 8+60.5/?i5+44.0/?i6+177.1/?3+49.5/? 6 +17.6/?ii . (1) 

Similarly, the required forces can be obtained for other load conditions using the analysis load 
combination in table 1 for the above-mentioned cases. For all combinations, the critical loads were 
found to be lift off and high-Q load conditions. 1 Figure 1 and table 1 give details of these conditions. As 
will be shown later, this method is intended to be used as one of the alternate approaches to simulate the 
load conditions for finite element analysis and the corresponding principal stresses without actually 
rerunning the finite model of the ETA ring. 

Using the standard techniques of theory of elasticity, the stresses were calculated determinis- 
tically for each ring component at critical sections identified by the NASTRAN output postprocessing. 

Net (minimum) section properties were used in these analyses based on the worst-case dimen- 
sions. This is the traditional minimum resistance-maximum load approach of deterministic analysis. The 
corresponding margins of safety (MS) are also calculated. 
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It is to be noted that all these calculations corresponding to the deterministic model are per- 
formed based on the assumption that the material and the load properties are also deterministic. Since in 
reality all variables have certain variation, small or large depending on the actual case, these should be 
considered in the analysis. Hence probabilistic analysis is performed to study the actual safety reserve. 


IV. BRIEF LITERATURE REVIEW 


There has been a significant amount of work done in the past two decades in the area of relia- 
bility. Notable of these are references 3 to 9 in terms of the basic development of reliability concepts. In 
addition, there have been several applications by various authors. 10-17 This report does not deal in detail 
with the work done by the above-referenced authors, but the work done for this report is a continuation 
along the same lines but applied to shuttle structure and specifically to the ETA ring. It should also be 
pointed out at this stage that the literature review quoted here relates to the work done for the present 
report and is by no means an exhaustive review of all the work done in the reliability area. 


V. RELIABILITY ANALYSIS 


To perform a reliability analysis, the formulation of limit state function g is essential. The corre- 
sponding safety index (/?) which is defined in subsequent paragraphs can then be calculated. 

In this study, the reliability index (/J) for a ring connecting the SRB to the ET for the limit state 
of stress is to be calculated. ft in turn depends on the mean value standard deviation of the random vari- 
ables connected with the material, geometric, and load properties. Specifically, /? can be defined as 


where g is the limit state function, g is the mean value, and a g is the standard deviation of g. g in turn 
can be expressed as 


g = R-S , (3) 

where R is the resistance and S is the load function for the structure under consideration. As can be seen 
in figure 8, p can be considered as the measure of the number of standard deviations that the mean value 
of the limit state function is from the failure surface. Figure 9 shows the demarcation of failure and 
survival states. 14 If R and S both are normal or lognormal, p and the corresponding probability of failure 
are related as 


or simply 


P = 0~\\-P f ) , 
Pf = 0(-p) , 


m 


m 


where <T>is the standard cumulative normal distribution function, and is the inverse of the standard 
cumulative normal distribution function. Figure 10 shows the actual p/ when R and S are both random 
variables. The reliability R can then be expressed as 
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Figure 10. Probability of failure measure. 

R = \-p f . (5) 

This implies that an increase in /} leads to an overall increase in the reliability of the structure. It should 
be noted that equations (4a) and (4 b) are not exact when the resistance and the load parameters are not 
normally distributed and/or when the limit state function is not linear in its variables. 

The problem would be simple if there exists a closed-form expression for the resistance variable, 
load variable, and the limit state function. Even in those cases, the exact evaluation of the probability of 
failure integral is almost impossible for non-normal cases. The general probability of failure function is 
given: 

P f = /•/ f x (x l ,x 2 ,.~,x n )dx 1 dx 2 ...dx n , (6) 

in which f x is the joint probability density function for x\,x 2 , .... x„ and the integration is performed over 
the region where g < 0. The above expression is based on a limit state function of g(xi, x 2 , ..., x„) = 0. In 
such cases, Monte Carlo simulation techniques are used to evaluate p/, R, or /?. 

In the present study, as the structure is complicated there are no closed-form relations between 
the output variable of stress and the corresponding input variables. Hence, probabilistic study is even 
more complicated in this case. The deterministic procedure used for the calculation of stresses is the 
well-known finite element method (FEM). The corresponding probabilistic study, which will use Monte 
Carlo techniques, will involve an FEM method in conjunction with the properties of the probability dis- 
tributions of all variables involved in the calculation of the stresses at critical points. 
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Two approaches were considered in this stage. The first approach was to treat the material and 
load properties as random variables while treating the geometric properties as deterministic. In this case, 
actual finite element analysis is to be performed for simulated load and material properties to give prin- 
cipal stresses using Monte Carlo simulation (discussed in subsequent paragraphs). The geometric prop- 
erties were to be considered as deterministic to avoid remodeling of this complex mesh, realizing that 
this is a shortcoming of this approach. The second approach is to treat the material, geometric, and load 
properties as random variables but utilizing the results of the deterministic finite element analysis. This 
is done by considering the variation of cross-sectional parameters, variation of grid point forces, mem- 
brane forces, and bending moments, etc. (output of deterministic finite element analysis), and variation 
of allowable stress of steel into the model and using Monte Carlo simulation to get the simulated 
principal stresses. It is to be noted that in this approach the load and material properties are not explicit 
random variables, but implicitly they are treated as random variables by treating the grid point forces 
and membrane forces, etc., as random variables. It is this second approach that is mainly used in this 
report. 


The procedure is explained with respect to the limit state of stress. For a general case, the prob- 
ability of failure for the limit state of stress can be expressed as: 

Pf = F(°act > Call) » (7 a) 

or 

Pf = P(<*ALL < <7aCt) • {lb) 

This is like pf= p(R < S) as traditionally used, where R is the combined resistance function, and 5 is the 
combined load function. The word combined is used here because <r, whether it is actual or allowable, is 
based on the material, geometric, cross-sectional, and load properties, and hence, Call and Cact are not 
clear cut resistance (where, in general, resistance function would involve only material, geometric, and 
cross-sectional properties) and load (would involve only load parameters). Cact can be expressed as: 

^ ACT — g 1 '^g2 > •••'^gn'^cX '^c2'"“>^cn'^L\ '^L2 * i) » (8) 


where 


x m \ = first variable representing material property 
x g \ = first variable representing geometric property 
x c \ = first variable representing cross-sectional property 

xli = first variable representing load properties directly or the output parameter of the finite 
element analyses (like the membrane force, grid point force, etc.). 

In the present study, <t A ct is the principal stress o P . The design nominal (deterministic) values of these 
variables were obtained from the structures materials/load group 1 at MSFC and from the available litera- 
ture. Since a probabilistic study is to be performed with the intent of final calculation of reliability, the 
variation of these parameters is required. To be specific, the statistical properties (like the mean value 
and coefficient of variation) of the distribution and the type of the distribution followed by these vari- 
ables is required. The relevant expressions are given as: 
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(9) 


X = BX , 

v.-pM 1 2 ]". (10) 

v B =[v?+vi+...v, x \ m , ( 11 ) 

where X and V x are the estimates of the mean value and the coefficient of variation based on the avail- 
able data; X and V x are the “true” mean and coefficient of variation of the variable; and B and V B are the 
mean value and standard deviation of bias factor. If the model is unbiased, B = 1 and Vg = 0. If the data 
are limited, jc and V x can be calculated from, 


x = 


x { +x 2 


( 12 ) 


_ *2— *1 
v x ~ 1 » 

J (x 2 +x ,) 


(13) 


where it is assumed that 95 percent of the values lie between x\ and x 2 . x\ and x 2 can be considered as 
the practical extreme limits. 


Similarly, (Tall can be expressed as. 


G ALL — fl(x m 1 yX m2 , ...fXfwfXgi ,Xg 2 ,. 


•’XgntX c j ,X C 2, •••>X cn ,X B i , X B 2 » ""*X \ n ) 


(14) 


which is a different function using the allowable values of the input variables. Alternatively, (Tall may 
be deterministic if it is a value taken from the governing codes for that structure. In any case, there is no 
closed-form expression for (Tact or (Tall- The computer will take the random values of the input vari- 
ables, go through finite element code (or use the results with a certain coefficient of variation), and using 
Monte Carlo techniques will calculate the stresses at different critical points. Once finite element analy- 
sis of these hypothetical rings is done (or the variation of the output parameters is considered) in con- 
junction with Monte Carlo techniques, a statistical analysis is done to calculate p/ and the safety index 
along with the reliability level using equations discussed above. Similar reliability analysis can be done 
for other limit states of the structure. In the present study, (Tall is (<t«)all the ultimate allowable or (cty) 
the yield allowable stress depending on the actual case being considered. 

The detailed procedure of Monte Carlo simulation in conjunction with the deterministic finite 
element analysis is given below: 18 

1 . Generate a random number r, between 0 and 1 using Rand Corporation procedure or stand- 
ard subroutines in the computer corresponding to a uniform distribution for any input variable dealing 
with cross section, load/force parameter and also for allowable stress. Note at this stage all the necessary 
random numbers for all input variables are generated either dealing with both actual stress (Tact or 0p 
and also allowable stress for (Tall- 


2. For the actual random value to be generated for any variable corresponding to a given distri- 
bution, the statistical properties of that variable are needed, like x and V x . Calculate x and V x from the 
available experimental data and then calculate x and V x from equations (9) to (1 1) using the proper bias 
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factor. In the absence of sufficient data for any of the input random variable x, calculate x and V x from 
equations (12) and (13). 

3. Since all the variables are assumed to be normal, the transformation of random numbers 
from uniform to normal is performed using. 


Ci = 0-Hn) , (15) 

where C,- is a standardized random number such that it is normally distributed with mean = 0 and vari- 
ance = 1. 0 l is the inverse of the standard cumulative normal distribution function. Note that standard 
subroutines are available to perform such an operation. IT is also to be noted that equation (15) uses the 
r, values generated in step 1. 

4. Knowing C, values from equation (15) and x and V x from step 2, a random value of x can 
then be generated using 


Xi = x+s x Ci , (16) 

where s x is the standard deviation of * given by s x = xV x . 

5. Generate random values for all independent variables (as mentioned in step 1) following 
steps 1 to 4 and then using Monte Carlo simulation and the equations of theory or elasticity to calculate 
principal stresses for all the critical elements of these hypothetical rings subjected to hypothetical 
loads/forces. It is to be noted that the word hypothetical is not used to indicate that these rings and the 
loads are all fictitious but instead it indicates that the deterministic ring problem is now studied by per- 
turbing its material, geometric, and load properties consistent with the actual structural behavior. Since it 
is now widely accepted that no variable is truly deterministic, this study attains significant importance. 

6. Having obtained a set of random values for each of the dependent variables like Gall 
Oact» perform statistical analysis of the required parameters. This could be the mean value, standard 
deviation, and coefficient of variation of all the response parameters. 

7. Calculate the probability of failure of an event corresponding to equation (7 b) from 


= T ■ (H) 

where rt\ is the number of times that an event has occurred in a sample size of n. For example, to calcu- 
late pf corresponding to equations (la) and (lb), once a random number has been generated for Call or 
<T P and Cact. it is checked if <t A ll < Oact- If it is true, then it corresponds to a failure event. Such occur- 
rences are totaled for the sample size, and pf is then calculated from equation (17). 

8. Knowing pf, one can calculate reliability R from equation (5). 

9. The safety index /3 can be calculated from equation (4 a) if the distribution of all the input 
variables is normal. Otherwise (3 can be calculated using equation (2) where g and <y g are the results 
obtained from step 6. However, this approach is not used in this study as all variables are assumed 
normal. 
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10. Steps 6 to 9 can be repeated for all the limit states of the ring problem. In the present study, 
only one limit state (stress) is considered. 

11. The effect of coefficient of variation of some of the critical input parameters on reliability of 
the ring (with respect to various limit states) will be studied as part of the overall behavioral study of the 
ring. 

If resistance and load functions are approximated as normal, then the safety index (/J) has a 
closed-form expression given by. 


/?-S (n-1) , 

\ln 2 Vl+Vf ' 


(18) 


where n = RJS , V R = a JR, and Vs = aJS. For the notations of this report, R is <t A ll and S is the ctact or 
Op. Similarly for resistance and load functions being lognormal, p is given by, 

' n(n) - . (19) 

•Jvj+vi 

For the variables not conforming to the above standard cases, p can be obtained numerically. As 
explained earlier, R and S in this report correspond to <7all and <t A ct. respectively. 


VI. FORMULATION OF THE PROBLEM 


As mentioned previously in this report, the components of the ETA ring that are critical were 
tunnel splices and H-fittings. In each of these components, there are few sections that were found critical 
based on deterministic analysis. 1 For all these cases, the reliability levels with respect to the limit state of 
stress are calculated. The final output of the probabilistic analysis is principal stress CTp. Hence, the prob- 
ability of failure corresponding to equation (lb) can be written as: 

Pfy = P((Gy) ALL < a p) * 


similarly for ultimate allowable stress 


/Vu = P((<Oall<<V • (21 

Knowing pfy and pf u , the corresponding reliability values R y and R u can be easily calculated from the 
generic equation (5) given earlier. The safety indices p can be calculated from equation (18) as the vari- 
ables have been assumed to be normal. 


At this point, it is important to note that in a deterministic analysis the MS is calculated using the 
equation given below: 1 


(A/S) yield - 


( g y) ALLOW 

1-1 (^p) ACTUAL 


( 22 ) 
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for ultimate allowable stress: 


(M5) ult - 


(g Fallow 

1-4 (O’ p) ACTUAL 


(23) 


The factors 1.1 and 1.4 that are used in equations (22) and (23) represent the safety factor in yield and 
ultimate allowable stress. Since safety factor is nothing but an ignorance factor and was introduced to 
compensate for all uncertainties, it should be discarded in reliability analysis as is traditionally done. 
That raises an important question with regard to equations (20) and (21). If in an analysis the data show 
that the basic variability in (<*>,) all and (Ou)all is the same, that would indicate that pfu is smaller than 
Pfy and hence is not as critical as the yield case. This is opposite to the general notion in deterministic 
analysis wherein the ultimate case is more critical than yield due to different safety factors. This should 
also be considered along with the fact that taking sf= s/S (where s is the allowable stress and 5 is the 
actual stress). It can be proved that 


considering both as random variables. In cases where either as or a s , the standard deviation in allowable 
and actual stress, is high, V s f, the coefficient of variation could be high. In some cases, this may have to 
be included in the model as an additional uncertainty in addition to the usual uncertainty in material 
geometric and load properties making sure it is not considered twice. On the whole, the authors of this 
report feel that so far as the reliability analysis is concerned it is the ultimate allowable stress (limiting 
stress) that is important. Provided that data are collected properly and proper statistical tools are used, 
the uncertainty in material, geometric, and load properties should take care of the safety factor auto- 
matically, and the reliability analysis should give a reasonable idea of the safety reserve in the structure. 
The results in this report are given using limit state stress of (<t«)all • 

For a ready comparison of probabilistic margin of safety (PMS), values are calculated for com- 
parison with the corresponding deterministic MS. This is given as 


PMS = 


(O' u) ALLT p) MAX 


J OU 


(24) 


where. 


(o p) max = Op+pS^p > (25) 

here, 

(0m)all = expected value of allowable stress 
S m = standard deviation of allowable stress 
Sop = standard deviation of principal stress 

a p = mean value of principal stress 

P = safety index of the structure with respect to stress limit state. 
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The reliability analysis of the above mentioned critical components will now be discussed in 

detail. 


The results for the basic cases of critical sections (cases A to E) are obtained for a sample of size 
of 5 000 random numbers using the variance reduction techniques 4 and for the following basic values of 
input parameters: V c = 0.01, V f = 0.06, V AS = 0.022, C MNC = U5, C MNF = 1.0, and C MNS = 1.055, where 
V stands for coefficient of variation, and C stands for coefficient of mean/normal. The subscripts C, F, 
and A s for V stands for cross-sectional parameters, forces (grid point forces, membrane forces, and 
bending moments), and allowable stress, respectively. The subscripts MNC, MNF, and MNS for C stand 
for mean to nominal for cross-sectional parameters, mean to nominal for forces (grid point, membrane, 
and bending moments), and mean to nominal for allowable stress parameters, respectively. The effects 
of other V and C values are studied under sensitivity analysis. The V c values are actually lower/higher 
than 0.01 for some cross-sectional parameters based on the general data in the literature. 12 16 20 To avoid 
any ambiguity, the probabilistic data for input parameters are tabulated for cases A to E when those 
cases are dealt with in the report individually. 


A. 90° Forward Tunnel Splice Plate (Minimum Section) 


This corresponds to USBI drawing 10170-0367 from reference 1. The general expression for 
principal stress is given by, 1 


<J x +0 y 

a p\> a pl~ 2 




+ X 




<T = MAX((T pl ,G p2 ) , 


(26) 

(27) 


which is considered for reliability calculation. 


x = — 

l *y A ’ 


Fy Mf, 

a y~A y + i xc 


_ F M z C y M£, 

0 *-A x + I ' + ly 


(28) 

(29) 

(30) 


Oi, <y y , and Xxy represent normal stresses in x and y direction and shearing stress, respectively. F# F y , Fxy, 
A/*, M y , and M z are the element membrane forces and bending moments. A x , A y , I xc , I y , and I z are cross- 
sectional properties. These parameters can be obtained knowing the critical cross section. The critical 
cross section for this is shown in figures 1 1 and 1 2. 

It is to be noted that all the values shown are nominal/design values. Table 2 shows the list of all 
the basic input variables used in the probabilistic analysis including that for allowable stress of steel. The 
information regarding x/xn for any variable and the coefficient of variation is obtained from discussions 
with the relevant groups at MSFC and based on existing data in the literature. 14 20 22 The results are 
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Figure 11. Cross-sectional parameters and grid point forces for critical section of case A (from ref. 1). 



Figure 12. Cross-sectional parameters for critical section of case A (ref. 1) (view AA of fig. 1 1). 
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Table 2. Case A: Probabilistic input parameters. 



obtained in terms of probability of failure using Monte Carlo simulation discussed earlier in this report. 
The corresponding reliability is then calculated from equation (5). The safety index p and PSM are cal- 
culated from equations (18) and (24), respectively. The results are shown below: 

a p = 91M Sop = 9.39 189.9 Sou = 4.17 

p fu = 0.15xl0- 14 R u = 0.9(H) P = 8.93 PMS = 1.89. 

This shows a high PMS compared to the corresponding deterministic MS of -0.03. It is to be noted that 
the principal stress has a high coefficient of variation. Figure 13 shows the Warner diagram for case A. 

A detailed discussion appears later. 
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DENSITY FUNCTIONS FOR SIGP1 AND SIGMUA 

SGUM= 1 74,SGUMX=204.SGP 1 M=7 1 ,SIGP 1 MX= 1 38 



Figure 13. Warner diagram for case A. 







B. 90° Aft Tunnel Splice (Small Section) 

This corresponds to USBI drawing 10170-0372 from reference 1. The critical section is shown 
in figure 14. The input data are shown in table 3. The results are shown below. 

131.39 S(jp = 14.98 a u = 189.9 5^ = 4.17 

Pfu = 8.03xl(H R u = 0.99919 p u = 3.7604 PMS = 0.5143. 


Again a good reliability level is seen in the results. The PMS of 0.5143 is much higher than the corre- 
sponding deterministic MS of -0.24. The corresponding probability of failure diagram is shown in 
figure 15. 
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b) CROSS SECTIONAL PROPERTIES AND 
GRID POINT FORCES 

Figure 14. Cross-sectional properties and grid point forces for critical section of case B (from ref. 1). 
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Table 3. Case B: Probabilistic input parameters. 


Parameters 

Nominal 

Mean/Nom 

Est. Cov. 

Bias 

Cov. Bias 

D1 

2.250 

1.15 

0.0042 

1.01 

0.05 

D2 

0.590 

1.15 

0.0042 

1.01 

0.05 

D3 

0.600 

1.15 

0.0042 

1.01 

0.05 

D4 

0.300 

1.15 

0.0042 

1.01 

0.05 

D5 

1.160 

1.15 

0.0042 

1.01 

0.05 

T21 

10,967.000 

1.0 

0.06 

0.95 

0.05 

T221 

6,243.000 

1.0 

0.06 

0.95 

0.05 

T222 

10,236.000 

1.0 

0.06 

0.95 

0.05 

T231 

6,681.000 

1.0 

0.06 

0.95 

0.05 

Fxl 

29,176.000 

1.0 

0.06 

0.95 

0.05 

Fx2 

28,134.000 

1.0 

0.06 

0.95 

0.05 

Mxl 

485.000 

1.0 

0.06 

0.95 

0.05 

Mx2 

453.000 

1.0 

0.06 

0.95 

0.05 

Fy 

-917.000 

1.0 

0.06 

0.95 

0.05 

Myl 

41.200 

1.0 

0.06 

0.95 

0.05 

Fxyl 

2,117.000 

1.0 

0.06 

0.95 

0.05 

Fxy2 

1,675.000 

1.0 

0.06 

0.95 

0.05 

FSY 

1.100 

1.0 

0.06 

1.00 

0.00 

FSU 

1.400 

1.0 

0.06 

1.00 

0.00 

SIGy 

163.000 

1.055 

0.022 

1.00 

0.00 

SIGu 

180.000 

1.055 

0.022 

1.00 

0.00 
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DENSITY FUNCTIONS FOR SIGP1 AND SIGMUA 

SGUM= 1 74,SGUMX=203.SGP1 M=96,SIGP1 MX= 1 88 



SIGP1 .SIGMUA 

□ SIGMUA 4- SIGP1 

Figure 15. Warner diagram for case B 




C. 90° Forward Tunnel Splice (Critical Section) 

This corresponds to USBI drawing 10170-0367 of reference 1. The critical sections are shown in 
figure 16. The probabilistic input parameters used for reliability analysis are shown in table 4. The 
results are shown below. 

e p = 118.27 Sop = 21.39 a u = 189.9 S^^AAl 

p fu = 0.0036 R u = 0.9964 p u = 3.285 PMS = 0.317. 


The corresponding deterministic MS is -0.22. This shows that probabilistic analysis has predicted a 
higher safety reserve than was calculated in the deterministic analysis. The corresponding Warner dia- 
gram is shown in figure 17. 




b) CROSS SECTIONAL PROPERTIES AND 
GRID POINT FORCES 

Figure 16. Cross-sectional properties and grid point forces for critical section of case C (from ref. 1). 
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Table 4. Case C: Probabilistic input parameters. 


Parameters 

Nominal 

Mean/Nom 

Est. Cov. 

Bias 

Cov. Bias 

T 

0.235 

1.15 

0.0400 

1.01 

0.05 

D1 

0.285 

1.15 

0.0042 

1.01 

0.05 

D2 

0.510 

1.15 

0.0042 

1.01 

0.05 

D3 

2.134 

1.15 

0.0042 

1.01 

0.05 

D4 

0.202 

1.15 

0.040 

1.01 

0.05 

D5 

0.769 

1.15 

0.040 

1.01 

0.05 

D6 

4.030 

1.15 

0.0042 

1.01 

0.05 

D7 

3.210 

1.15 

0.0042 

1.01 

0.05 

D8 

1.960 

1.15 

0.0042 

1.01 

0.05 

D9 

1.800 

1.15 

0.0042 

1.01 

0.05 

DIO 

0.820 

1.15 

0.0042 

1.01 

0.05 

Dll 

1.250 

1.15 

0.0042 

1.01 

0.05 

D12 

1.350 

1.15 

0.0042 

1.01 

0.05 

D13 

0.740 

1.15 

0.0042 

1.01 

0.05 

D14 

1.660 

1.15 

0.0042 

1.01 

0.05 

D15 

1.040 

1.15 

0.0042 

1.01 

0.05 

Wm 

2.450 

1.15 

0.0042 

1.01 

0.05 

T21 

7,498.000 

1.0 

0.060 

0.95 

0.05 

T221 

10,470.000 

1.0 

0.060 

0.95 

0.05 

T222 

10,749.000 

1.0 

0.060 

0.95 

0.05 

T231 

15,603.000 

1.0 

0.060 

0.95 

0.05 

T232 

13,299.000 

1.0 

0.060 

0.95 

0.05 

T241 

11,100.000 

1.0 

0.060 

0.95 

0.05 

T242 

5,717.000 

1.0 

0.060 

0.95 

0.05 

Fy 

1,077.000 

1.0 

0.060 

0.95 

0.05 

Mxl 

24.100 

1.0 

0.060 

0.95 

0.05 

Mx2 

75.600 

1.0 

0.060 

0.95 

0.05 

Mx3 

-8.900 

1.0 

0.060 

0.95 

0.05 

Mx4 

-6.700 

1.0 

0.060 

0.95 

0.05 

Myl 

237.000 

1.0 

0.060 

0.95 

0.05 

Fxyl 

-1,005.000 

1.0 

0.060 

0.95 

0.05 

Fxy2 

-2,539.000 

1.0 

0.060 

0.95 

0.05 

Fxy3 

-399.200 

1.0 

0.060 

0.95 

0.05 

Fxy4 

402.000 

1.0 

0.060 

0.95 

0.05 

Fx 

5,476.000 

1.0 

0.060 

0.95 

0.05 

Fsy 

1.100 

1.0 

0.01 

1.00 

0.00 

Fsu 

1.400 

1.0 

0.01 

1.00 

0.00 

SIGy 

163.000 

1.055 

0.022 

1.00 

0.00 

SIGu 

180.000 

1.055 

0.022 

1.00 

0.00 
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DENSITY FUNCTIONS FOR SIGP1 AND SIGMUA 

SGUM= 1 74,SGUMX=203,SGP1 M=65,SIGP1 MX=21 6 



SIGP1, SIGMUA 

□ SIGMUA + SIGP1 

Figure 17. Warner diagram for case C 








D. H-Fitting (Plates) 


This corresponds to USBI drawing 10170-0011 from reference 1. The critical section and the 
associated cross-section parameters are shown in figures 18 and 19. The probabilistic input parameters 
are shown in table 5. The results of reliability are shown below. 

a p = 72.396 5^ = 6.503 d u = 189.9 5^ = 4.177 

p fu = (<10- 14 ) R u = 1.09 P u = 15.201 PMS = 4.461 . 

It is to be noted that the PMS is quite high. The corresponding deterministic MS was +0.52. The corre- 
sponding Warner diagram is shown in figure 20. 


G. H-Fitting (Lugs) 

This corresponds to USBI drawing 10170-011 of reference 1. The critical section and the asso- 
ciated cross-sectional parameters are shown in figures 21 and 22. The probabilistic input parameters are 
shown in table 6. The results are shown below. 

o p = 116.15 Sap= 14.39 a„= 189.9 So* = 4.17 

Pfu = 0.37x10-0 R u = 0.9(6) Pu = 4.920 PMS = 0.699. 

Once again it has been shown that PMS (0.699) is much higher than the deterministic MS of -0.21. The 
corresponding Warner diagram is shown in figure 23. 


VII. SENSITIVITY ANALYSIS 


As mentioned above, cases A to E were studied for some basic core parameters as discussed in 
section VI. Since no actual experimental data were available for the input parameters and the data used 
were based on the available data in the literature and the discussions with the load and stress group at 
MSFC, it was found prudent to study the effect of variation of the input parameters on the reliability of 
the ETA ring for the stress limit state for most critical components at critical sections. Based on the 
results obtained for cases A to E and looking at the corresponding Warner diagrams, it is seen that the 
overlap between the actual principal stress and the allowable ultimate stress is higher for cases B and E 
compared to cases A, C, and D; while the PMS is consistently higher than the corresponding MS for 
deterministic analysis. Hence, the sensitivity analysis study is conducted for case B (tunnel splice plate) 
and case E (H-fitting lugs). 

This study is done by studying the effect of variation of coefficient of cross-sectional parameters 
(V c ), force parameters (V/) like membrane forces, bending moments, etc., allowable stress parameters 
(V A5 ) on the reliability of the ETA ring, taken one at a time. The basic values considered are V c = 0.01, 
V f = 0.06, and V AS = 0.022. These results are shown in figures 24 to 27 for case B and in figures 28 to 31 
for case E. In addition to this, the effect of X/xn denoted as c mn (coefficient of mean to nominal), on 
reliability is studied for various cross-sectional parameters (C mnc ), force parameters {C mn f), and 


29 




sec 









Table 5. Case D: Probabilistic input parameters. 


Parameters 

Nominal 

Mean/Nom 

Est. Cov. 

Bias 

Cov. Bi 

T1 

0.430 

1.15 

0.0400 

1.01 

0.05 

D1 

0.515 

1.15 

0.0042 

1.01 

0.05 

D2 

0.971 

1.15 

0.0042 

1.01 

0.05 

D3 

0.421 

1.15 

0.0042 

1.01 

0.05 

D4 

1.366 

1.15 

0.0042 

1.01 

0.05 

D5 

0.964 

1.15 

0.0042 

1.01 

0.05 

D6 

0.325 

1.15 

0.0042 

1.01 

0.05 

D7 

6.940 

1.15 

0.0042 

1.01 

0.05 

D8 

5.820 

1.15 

0.0042 

1.01 

0.05 

D9 

4.470 

1.15 

0.0042 

1.01 

0.05 

DIO 

3.395 

1.15 

0.0042 

1.01 

0.05 

Dll 

2.123 

1.15 

0.0042 

1.01 

0.05 

D12 

1.328 

1.15 

0.0042 

1.01 

0.05 

D13 

0.163 

1.15 

0.0042 

1.01 

0.05 

YBN 

3.345 

1.15 

0.0042 

1.01 

0.05 

D14 

6.395 

1.15 

0.0042 

1.01 

0.05 

D15 

5.145 

1.15 

0.0042 

1.01 

0.05 

D16 

3.795 

1.15 

0.0042 

1.01 

0.05 

D17 

2.445 

1.15 

0.0042 

1.01 

0.05 

D18 

1.525 

1.15 

0.0042 

1.01 

0.05 

D19 

0.585 

1.15 

0.0042 

1.01 

0.05 

D20 

0.940 

1.15 

0.0042 

1.01 

0.05 

D21 

0.920 

1.15 

0.0042 

1.01 

0.05 

D22 

1.350 

1.15 

0.0042 

1.01 

0.05 

D23 

1.250 

1.15 

0.0042 

1.01 

0.05 

Wm 

2.150 

1.15 

0.0042 

1.01 

0.05 

D24 

7.200 

1.15 

0.0042 

1.01 

0.05 

Wmin 

1.740 

1.15 

0.0042 

1.01 

0.05 

T21 

10,252.000 

1.0 

0.060 

0.95 

0.05 

T221 

17,205.000 

1.0 

0.060 

0.95 

0.05 

T222 

14,139.000 

1.0 

0.060 

0.95 

0.05 

T231 

15,936.000 

1.0 

0.060 

0.95 

0.05 

T232 

15,530.000 

1.0 

0.060 

0.95 

0.05 
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Parameters 


T241 

T242 

T251 

T252 

T261 

Fyl 

Fxyl 

Fxy2 

Fxy3 

Fxy4 

Fxy5 


Mean/Nom 


Est. Cov. 


Cov. Bias 


8.955.000 

11.708.000 

3.473.000 

13.460.000 

4.146.000 

9.298.000 
-3,531.000 
-6,824.000 

-11,864.000 

-11,186.000 

-9,245.000 

-82.100 


SIGy 

SIGu 


-230.000 

-625.000 

-440.000 

0.250 

1.100 

1.400 

150.000 

180.000 


1.055 

1.055 
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Figure 20. Warner diagram for case D. 











a) CROSS SECTIONAL PROPERTIES 



b) CROSS SECTIONAL PROPERTIES AND 
GRID POINT FORCES 

Figure 22. Cross-sectional properties and grid point forces for critical section of case E (from ref. 1 ). 
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Table 6. Case E: Probabilistic input parameters. 


Parameters 

Nominal 

Mean/Nom 

Est. Cov. 

Bias 

Cov. Bias 

TH1 

1.440 

1.15 

0.0400 

1.01 

0.05 

D1 

0.760 

1.15 

0.0042 

1.01 

0.05 

D2 

0.770 

1.15 

0.0042 

1.01 

0.05 

D3 

0.990 

1.15 

0.0042 

1.01 

0.05 

D4 

1.499 

1.15 

0.0042 

1.01 

0.05 

Til 

-41.300 

1.0 

0.06 

0.95 

0.05 

T21 

30.060 

1.0 

0.06 

0.95 

0.05 

T121 

9.930 

1.0 

0.06 

0.95 

0.05 

T122 

- 6.720 

1.0 

0.06 

0.95 

0.05 

T221 

58.430 

1.0 

0.06 

0.95 

0.05 

T222 

53.520 

1.0 

0.06 

0.95 

0.05 

T13 

0.780 

1.0 

0.06 

0.95 

0.05 

T231 

72.500 

1.0 

0.06 

0.95 

0.05 

Mxl 

4.680 

1.0 

0.06 

0.95 

0.05 

Mx2 

0.140 

1.0 

0.06 

0.95 

0.05 

Myl 

0.260 

1.0 

0.06 

0.95 

0.05 

My2 

3.920 

1.0 

0.06 

0.95 

0.05 

P 

214.510 

1.0 

0.06 

0.95 

0.05 

FSY 

1.100 

1.0 

0.01 

1.00 

0.00 

FSU 

1.400 

1.0 

0.01 

1.00 

0.00 

SIGy 

150.000 

1.055 

0.022 

1.00 

0.00 

SIGu 

180.000 

1.055 

0.022 

1.00 

0.00 
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DENSITY FUNCTIONS FOR SIGP1 AND SIGMUA 

SGUM=> 1 74,SGUMX=20 J.SGP 1 M«=80,SIGP 1 MX=* 1 76 



SIGP1, SIGMUA 

□ SIGMUA + SIGP1 

Figure 23. Warner diagram for case E. 







CASE B 



(ny) A±ni8vn3y 
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Figure 27. Combined effects of coefficient of variation of input parameters (V c , Vf, and V^y). 





CASE E 



(ny) Aimavraa 
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Figure 29. Effect of coefficient of variation of force parameters (Vp). 




CASE E 



Figure 31. Combined effect of coefficient of variation of input parameters. 




allowable stress parameters (C mns ). To be specific, C mnc stands for coefficient of mean to nominal for 
cross-sectional parameters, C mn f stands for coefficient of mean to nominal for force parameters, and Cmns 
stands for coefficient of variation of mean to nominal for stress parameters. The basic values considered 
are C mnc = 1.15, C mn f = 1.0, and C mns 1.055. These results are shown in figures 32 to 35 for case B and 
in figures 36 to 39 for case E. 

The effect of the sample size of random numbers on the reliability of the ETA ring is also studied 
by varying the number of random numbers from 5,000 to 20,000 and using varying reduction techniques 
as well as direct Monte Carlo simulation. The effect is found to be insignificant so far as the reliability 
values of the ETA ring are concerned. 


Vffl. DISCUSSION OF RESULTS AND CONCLUSIONS 


It can be seen from various cases discussed earlier in section VI that the coefficient of variation 
of the principal stress is greater than or equal to 0.10, and, in case D, it is as high as 0.18 even though the 
input coefficient of variation is not very high. This is because the expression for principal stress 
(equation (27)) is highly nonlinear and is a function of other dependent variables like normal and shear 
stresses. It should be noted that these stresses are again nonlinear functions of some dependent and inde- 
pendent variables. The cumulative effect of all this is reflected in the coefficient of variation of one of 
the limit state variables of principal stress. This in effect affects the reliability levels. 

Looking at figures 27 and 31, it can be seen that the V c has maximum effect on reliability for 
cases B and E. Similarly, looking at figures 35 and 39, it can be seen that Cmnc has maximum effect on 
reliability for cases B and E. Another interesting point to be noted from these figures is that the optimal 
value (the common intersection point) of C mn is approximately 1.1 for both cases B and E. 

On the whole, the technique of Monte Carlo simulation with and without variance reduction 
techniques has been used successfully for evaluation of reliability levels of the ETA ring for the stress 
limit state. The cross-sectional parameters are found to be most critical as their effect on reliability is 
found to be maximum. Reliability levels in general are found to be high and appear to be supportive of 
the deterministic results and with a higher safety reserve. 
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CASE B 
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Figure 33. Effect of coefficient of mean/nominal for force parameters (Cmnf)- 


CASE B 
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Figure 34. Effect of coefficient of mean/nominal for stress parameters (Cmns)- 



CASE B 



Figure 35. Combined effect of coefficient of mean/nominal parameters (C^vc. Q/jvf, and Cmns)- 
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